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Abstract: We explore a class of holographic superconductors built using non-abelian con- 
densates on probe branes in conformal and non-conformal backgrounds. These are shown 
to exhibit behaviour of the specific heat which resembles that of heavy fermion compounds 
in the superconducting phase. Instead of showing BCS-like exponential behaviour, the 
specific heat is polynomial in the temperature. It exhibits a jump at the critical tempera- 
ture, in agreement with real-world superconductors. We also analyse the behaviour of the 
energy gap and the AC and DC conductivities, and find that the systems can be either 
semi-conducting or metallic just above the critical temperature. 
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1. Introduction and summary 

An intriguing new application of the correspondence between string theory and gauge 
theory is the study of strongly coupled systems which exhibit superconductivity. Apart 
from suggesting an alternative point of view on non-BCS superconductors, these might also 
help to develop a better understanding of quantum critical points and strongly coupled 
superconductors for which so far a satisfactory theory is lacking [1]. Moreover, these 
constructions provide an interesting new arena in which to explore and test the string/gauge 
theory correspondence. We refer the reader to [2, 3] for more background information and 
references. 

In the original model of [1], a holographic superconductor was constructed by con- 
sidering an Abelian Higgs model in the background of an anti-de-Sitter black hole. A 
stable condensate of the scalar field can occur because of nontrivial coupling of the scalar 
to the gauge potential, in combination with nontrivial asymptotics [4] or because of the 
nontrivial near-horizon geometry which destabilises the asymptotically stable mode [1]. 
This system exhibits an energy gap which scales linearly with the critical temperature, 
and a frequency-dependent conductivity which shows some similarity with real-world su- 
perconductors. A full analysis involving back-reaction of the Abelian-Higgs fields on the 
gravitational background has also been carried out [5]. 
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A different way to obtain a non-trivial condensate in a holographic setup is to consider, 
instead, a non-abelian gauge field coupled to gravity. Again, because of the non-trivial 
asymptotics, it is possible to form condensates of the gauge field. For some time now, 
it has been known that such configurations exist even in systems at zero temperature 
(i.e. without black hole horizon), where they are related to vector meson condensates in 
the dual gauge theory [6]. An application to p-wave superconductors was developed in [4], 
based on earlier work in [7]. We will call these 'colourful superconductors' (not to be 
confused with colour superconductors in QCD). 

Many properties of these colourful superconductors have so far not been analysed, and 
it is far from clear which real-world superconductors they resemble most. Moreover, the 
analysis so far has focussed on conformal cases, since there is an expectation that ideas 
of the gauge/gravity correspondence are perhaps most useful in the context of quantum 
critical phenomena, which are often characterised by relativistic conformal symmetry. 

In the present paper we thus set out to investigate a variety of models which exhibit 
non-abelian condensates. We will focus mostly on D-brane constructions, both conformal 
and non-conformal, including the Sakai-Sugimoto model. These models obviously have 
an origin in a full-fledged string theory, however we will use them here without worrying 
about the validity of string theory approximations. Our attitude is pragmatic, in the sense 
that we are mainly interested in seeing how the models behave, with the eventual goal of 
improving the description of holographic superconductors and bringing them closer to real- 
world systems. For this reason, we will e.g. work with Yang-Mills action, without attempt 
to include DBI corrections to it, and we will also not consider gravitational back-reaction. 

After an introduction to the models under consideration, we will first focus on their 
thermodynamic properties, in particular the specific heat. For all systems which we anal- 
yse, we find evidence that in the superconducting phase, they resemble non-BCS, heavy 
fermion compounds. That is, their specific heat has a polynomial behaviour as a function 
of temperature, in contrast to the exponential one seen in BCS superconductors. At a more 
quantitative level, we find that c v ~ T n where n ~ 2.7 — 5.5 depending on the model. We 
also demonstrate the appearance of a jump in the specific heat at T c , a property common 
to both strongly and weakly coupled superconductors. Above T c , our approach only yields 
linear scaling for the specific heat at constant chemical potential for the D3/D7 system, 
while this is present also in the D4/D8 case in the DBI approximation [8]. The specific 
heat at constant density is, for all conformal cases, known to behave non-linearly [9]. 

We then continue in §|] with the analysis of electromagnetic properties. For most sys- 
tems under consideration we find behaviour above the critical temperature which resembles 
that of a semi-conductor, i.e. a decreasing resistivity with increasing temperature. In the 
superconducting phase, there is a small region in temperature in which a condensate has 
already formed but the resistivity is not yet entirely zero. We find that the gap scales 
linearly with the critical temperature, uj g oc T c , but the proportionality constant is out- 
side the BCS range and can also be substantially different from the ones found so far for 
conformal holographic superconductors. 

In two sections (§ |3.3| and § |4.3j ) we comment on the relation of our findings to real-world 
superconductors. 
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2. Colourful superconductors 



2.1 Non-abelian condensates 

Before we analyse a number of concrete holographic superconductors, let us first summarise 
the general setup and some universal aspects of the problem. We will consider generic 
diagonal metrics in which the time-time component of the line element is of the form 
—g(u) dt 2 , where g(u) has a zero at the location of the horizon. An example of this sort is 
the metric of a Dp-brane, given by 

dsl p = Q v (- f p (u)dt> + + (I) ^ (|fy + ^ 2 dnl-p) > (2-i) 

where i,j = 1 . . .p and f p (u) = 1 — (ut/u) 7 ~ p . For our numerical analysis later in the paper 
it is advantageous to introduce coordinates in which the coordinate distance between the 
asymptotic boundary and the horizon is finite. We will use r := L 2 /u. The Dp-brane 
metric then reads 

dsl p = (?) 2 (- f p (r)dt 2 + ^dsW) + (?) 2 (-^ + ? dnj_ p ) , (2.2) 

with f p (r) = 1 — (r/r^) 7_p . The horizon is now located at r = vt and the asymptotic 
boundary is at r = 0. In this coordinate system the dilaton equals 

T \ (p-7)(p-3)/4 

In addition to these Dp-brane systems we will also consider the AdS4 black hole background, 
with vanishing dilaton and metric given by 

d4ds 4 = (?) 2 (-fs(r)dt 2 + aydarW + |fy) . (2.4) 

For future reference, we note here that the temperature of the Dp-brane backgrounds, 
as measured by an asymptotic observer following the orbit of the Killing vector d/dt, is 
given by 

7 -^ L r. (2.5) 



47T L \VT J 

For the AdS4 black hole the temperature is T = 3/(47rry). 

If we now add a Dg-brane to the Dp brane background, the non-abelian gauge field on 
the Dg-brane can condense [6] ; a similar story holds true when we add a Yang-Mills action 
to the AdS4 background [7]. We will here restrict attention to an SU(2) gauge field. The 
non-abelian condensate will occur if the chemical potential [i, corresponding to the value of 
the Aq component on the boundary, is fixed to a large enough value. The system actually 
admits condensates of two different types [6, 7]. The one which generically turns out to 
have the largest energy of the two is of the form 

A = A { Q ] a 3 dt + A(r) {a l dx l + a 2 dx 2 ) . (2.6) 
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This solution breaks the SU(2) symmetry to a diagonal U(l). It is unstable [4], and 
expected to decay to a lower energy condensate. The lower energy configuration can be 
written in the form 1 

A = A { 3) a 3 dt + A { ^a l dx 3 . (2.7) 
Such a condensate is analogous to the one studied at zero temperature in [6], and we will 



return to this analogy in more detail in section 2.S. In the 2 + 1 dimensional case studied 
in [4], this condensate breaks both SU(2) and rotational invariance completely. 

We will restrict ourselves to an analysis of the Yang-Mills truncation of the effective 
action on the probe brane (comments on the limitations of this approximation will be made 
in § |3.2| ) . The Lagrangian which determines the formation of a condensate of the type ( |2.7| ) 
is given by 

L = -T q dr yf^e-* \g°°9 rr (d r A®f + g^{d r A^f + 4g™ g 33 (A^ A^ f] . (2.8) 
Jo 1 J 

This action leads to a coupled system of equations of motion, 



d r 

Or 



^ge-*g™rd r Af = A{^?Af y/^T*^ 



(2.9) 



These are to be supplemented with boundary conditions at infinity, which take the form 

A® = fi + 0(r) , = °( r ) ■ ( 2 - 10 ) 

and boundary conditions at the horizon, which follow from the requirement that A ' (ry) = 0. 
One typically finds that for values of the chemical potential below a critical one, A^ is 
identically zero and the condensate is purely abelian. For chemical potentials larger than 
[i = /x c , the A^ component turns on, making the condensate non-abelian. These solutions 
are parameterised by particular combinations of values of d r A^ and A% at the horizon. 
If one keeps increasing ju substantially beyond the critical value, new branches of solutions 
occur. These typically have a larger energy and hence are expected to be unstable. An 
example of the structure of these solutions is given in figure [l|. 

2.2 Overview of the models 

We will consider both three-dimensional and four-dimensional superconductors of the type 
described above. For the three-dimensional models will will focus on the AdS4 black hole 
and the non-conformal D2/D6 intersection. 2 Our main emphasis in the context of four- 
dimensional superconducting systems will be the one obtained from the Sakai-Sugimoto 
model [11]. However, in order to appreciate the variety of ways in which this model is a 



1 The use of the coordinate x 3 is a bit misleading, and should not be understood to imply that this 
condensate only occurs in models which have at least three space dimensions. 

2 The non-abelian condensate in the AdS4 black hole background was previously constructed in [7]. The 
D2/D6 system allows, to some limited extent, an analysis which includes back-reaction effects [10], and it 
would be interesting to investigate how these influence the results obtained in the present section. 
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Figure 1: The structure of the condensate solution space for the D3/D7 system, showing 
the chemical potential p as a function of the two parameters A' and B which determine 
the solution at the horizon. Clearly visible are the various branches of condensates. 

non-trivial superconductor, we will here also review some known facts about the D3/D7 
superconductor, and extend them with several new results. 

A Frobenius analysis of the equations near infinity reveals that the asymptotic be- 
haviour of A^ and A3 is 



AdS 4 : 


,(3) 

A o 


= M~ 


pr^ , 


A? 


= p! + p'r + ■■■ 


D2/D6 : 


,(3) 
^0 


= M- 


pr 2 + • ■ ■ , 




= H* + p'r 2 + • • 


D3/D7 : 


,(3) 
A 


= M~ 


pr 2 + --- , 


^3 


= p + p'r 2 + • • 


D4/D8 : 


A (3) 
^0 


= M~ 


pr 3/2 + ... 




= p' + p'r^ 2 + 



(2.11) 



In all of the cases we will be interested in solutions for which the chemical potential for 
A^ is turned on, while it is set to zero for A^; in other words we are looking for solutions 
which satisfy the boundary conditions p ^ and p! = 0. As usual we also require that, 
for regularity of the solution, at the horizon A (ry) = 0. The equations of motion then 
imply that the near-horizon behaviour of the fields is (for all cases) 

Af =A'{r-r T ) + --- , A® = B + • • • . (2.12) 

By performing a two-dimensional scan through the parameter space spanned by A' and B 
and selecting those values for which the solutions to (^1]) satisfy the boundary condition 
that p! = 0, we arrive at the condensate solutions. 

It is convenient to express these solutions in terms of a dimensionless condensate p, 



- 5 - 



0' //U A 2 




Figure 2: The dimensionless condensates p' for (left to right, top to bottom) the AdS4 
black hole, D2/D6, D3/D7 and D4/D8 systems, as a function of the dimensionless ra- 
tio T/fj, (or to be numerically precise, 5 ^ 2 //i). The dots indicate selected states for 
which further properties will be analysed in section ^. 



given for the various systems by 



AdS 4 : 


P = 


p/p 2 , 


D2/D6 : 


P = 


P/P 7/3 


D3/D7 : 


P = 


p/p 3 , 


D4/D8 : 


P = 


P/P 4 , 



(2.13) 



(similar expressions hold for p'). The advantage of using p and p' is that these only depend 
on one dimensionless quantity p/T, not on the product TL (see appendix [A| for more 
details). 

Plots of the dimensionless condensate p' versus the ratio T/ p are given in figure [|. We 
see that at fixed temperature the condensate is nonvanishing when p < p c , while p' = 
for p > p c . In other words, the system undergoes a second-order phase transition as the 
chemical potential is changed. By numerically fitting the curves we see that the critical 
coefficient a, which determines the scaling behaviour of the condensate in the vicinity of 
the critical point according to p ~ (T/p c — T/p) a , is a = 1/2, in agreement with Landau- 
Ginzburg theory. 
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Figure 3: Comparison of the energy of the two non-abelian condensates and the abelian 
one, for the D3/D7 system, as a function of the chemical potential. The "three-component 
condensate" (2.6) (middle curve) always has a larger energy than the "two-component 
condensate" (2.7) (lower curve), and both have, for fi > fi c w 2.0, a lower energy than 



the abelian one (upper curve). The other systems analysed in this paper exhibit similar 
behaviour. 



For the alternative condensate given in (p.6[), which is governed by the equations 



Or 



r 



8A'A, 



2 A ( 3 ) 







(2.14) 



the density as a function of the inverse chemical potential behaves similarly. However, 
the energy of this configuration, obtained by evaluating E = —L with L given by (|2.£| ) 
is consistently higher. A comparison is displayed in figure [3|. We will henceforth restrict 
ourselves to the lowest-energy condensate, which is the true ground state of the model. 

2.3 Comments on the zero temperature limit 

The strict zero temperature limit of holographic superconductors remains somewhat of a 
mystery. In principle, there are three different scenarios. In the first one, one takes a 
naive zero temperature limit of the black hole metric, and obtains a metric with a Poincare 
horizon. This is what happens in the D3/D7 system upon taking — ► 0. The resulting 
system suffers from the fact that a condensate now becomes singular, even when using the 
full DBI action. 

Another possibility is that the limit should be taken in such a way that an extremal 
horizon is obtained. Concrete examples of this type have not yet been found, but they would 
presumably avoid the triviality of the Poincare horizon case discussed above. An attempt 
to analyse the zero-temperature situation was made in [12] by constructing a solution 
interpolating between two different AdS4 metrics. However, because of the conformal 
symmetry in the infrared, the solution found there exhibits a power-law scaling of the 
conductivity for small frequency, and hence no gap [13]. 
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There remains on further option, which seems more natural in e.g. the Sakai-Sugimoto 
model. In contrast to e.g. the A0IS4 black holes used in many studies of holographic su- 
perconductors, the Sakai-Sugimoto model has a well-understood limit to zero temperature. 
For sufficiently low temperature, the system undergoes a first order phase transition [14]. 
One then obtains a D4/D8 system in which there is no horizon on the brane, and the 
D8-branes which were previously intersecting the horizon are now connected in a smooth 
way. This system admits non-abelian condensates (as analysed in [6]). The absence of a 
horizon automatically implies that there is no dissipation in the system. 



3. Thermodynamic properties 
3.1 Free energy and specific heat 

Let us now turn to an analysis of some of the thermodynamic properties of the super- 
conductors. We will first present the results for the holographic models and then we will 
compare these with real world superconductors in section |3.3| . 

There has been some confusion in the literature as to the parameters which should be 
kept fixed when computing the specific heat. From the point of view of an experimental 
setup, keeping the density p constant is the most natural. This type of computation has 
been done for the abelian condensate of the D3/D7 case, using a full DBI analysis, in [9]. 
Alternatively, one could consider keeping the chemical potential p constant, which is a 
natural thing to do from the point of view of some holographic setups. A computation 
of this type can be found in e.g. [8] . Depending on what is kept fixed, the results can be 
quite different in as far as the dependence on the temperature is concerned. We will here 
elucidate these differences and then continue to show how the results change when one 
considers a non-abelian condensate instead. 

The main technical ingredient in the computation of the temperature dependence of 
the free energy is the scaling symmetry discussed in section [A]. It allows us to study 
the temperature dependence by analysing the dependence on the chemical potential fi 
at fixed temperature. By writing the Euclidean action in dimensionless coordinates and 
dimensionless fields Aq and A3, we find that 

AdS 4 : S E (T) = (TLf x S E (ji/T) , 

D2/D6 : S E (T) = {TL)l x S E (ji/T) , 

(3.1) 

D3/D7 : S E (T) = (TLf x S E (fi/T) , 

D4/D8 : S E (T) = (TLf x S E (n/T) , 

In writing these factors we have taken into account that the integral over the Euclidean 
time circle produces a factor of 1/T. 

For the computation of the heat capacity we now use the free energy F = TS E together 
with 

ds dF 
c v = T— - with s = — -— - , (3.2) 
dT dT' K ' 
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Figure 4: Free energies and specific heat for the D3/D7 system, as a function of tem- 
perature, for fixed chemical potential \x (and hence with p = p(T)). The black curve 
denotes the abelian condensate with only Aq non-zero, while the red curve denotes the 
non-abelian condensate. 



where s is the entropy density, and we have suppressed volume factors. Using the form ( |3.1| ) 
it is easy to compute quantities like the free energy or the specific heat at fixed chemical 
potential fi. In order to compute at fixed density p, one has to express pt/T in terms of p 
using condensate curves like those of figure |2|. 

Figure |] shows curves for the free energy F(T) at constant pi\ the curves at constant p 
are qualitatively similar. After normalising F(T = 0) = 0, a double-logarithmic plot (not 
shown here) leads to practically straight curves, which suggests a fit of the free energy 
to a function of the form a + bT c (as opposed to something involving exponentials of the 
temperature). We will comment below on how this connects to real- world superconductors. 
From the fit of the free energy we then obtain c v by applying ([^) to the fit, and arrive at 
the following numerical estimates in the normal and superconducting phases 3 , 

normal phase superconducting 
C v (pb = ct.) c v (p = ct.) c v (p = ct.) 

AdS 4 : T° T 2 - 7 T 2 - 9 

D2/D6: T\ T 30 T 32 (3 ' 3) 

D3/D7: T T 38 T 39 

D4/D8: T 2 T 5 - 5 T 52 

As explained below, the results for the free energy in the abelian phase at constant p 
are unreliable in the Yang-Mills truncation, and we have therefore not indicated values 
for c v (p = const.) in the normal phase. The systematic error bars on the numerical expo- 
nents receive contributions from various stages of the construction and should thus only 
be used as indicators about the numbers of degrees of freedom. 

The only previous result about the specific heat of a non-abelian condensate is, as far as 
we are aware, given in the work of [15]. These authors claim that c v {p = const.) ~ T 3 for the 



3 Numerically it is more reliable to fit F rather than c v . The ansatz is consistent with the fact that c v 
should vanish as T — > 0, as required by the third law of thermodynamics. 
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Figure 5: Specific heat for the D4/D8, AdS4 and D2/D6 systems respectively; curves 
as in figure ||. The free energy plots of these systems are similar to figure |Ja, the only 
difference being the exponents of the curves. 



D3/D7 system (analysed using the non-abelian DBI action in two different prescriptions). 
Our computations for D3/D7 lead to a somewhat larger exponent: they also contain a 
factor T 3 , arising from the overall prefactor in ( |3.1| ), but we find an additional temperature 
dependence through S(p/T). This difference may be a consequence of the DBI corrections 
to the Yang- Mills action, but we did not check the results of [15]. 

The most important conclusion to draw from our analysis - more important than the 
actual numbers - is the fact that we find no evidence for exponential suppression of the 
specific heat in the superconducting phase. As we will see, this is in sharp contrast with 
BCS-like superconductors. The second important result is that the non-abelian systems 
exhibit a jump in the specific heat at the phase transition. We will discuss the relevance 



of these findings to real world superconductors in §3.2. 



3.2 Limitations of the Yang-Mills approximation 

The result for the specific heat at constant density for the abelian condensate (i.e. above 
T c ), which were obtained from our Yang-Mills analysis, differ substantially from the ones 
obtained using a DBI analysis [9]. However, we have also seen that computing the specific 
heat above T c at constant chemical potential, gives a much better agreement with the 
DBI computation [8]. The reason for this seems directly related to the behaviour of the 
corresponding two gauge field solutions in the zero-temperature limit, as the one at constant 
chemical potential admits a smooth T — > limit while the solution at fixed density does 
not. 

To see this explicitly, consider the abelian solution to the Yang-Mills equations of 
motion, for instance for the D4/D8 model. It reads 

d r Ao = prV 2 . (3.4) 

The chemical potential and the free energy are (making use of ( |2.5| ), i.e. T = 3/(4-7ry / rr)) 

rr T rr T 2/4 x ~ 3 

fi = / 8 r A = / drpr 1 / 2 = -p [ -ttT 



J 



3' V3 

(3.5) 



Jo Jo 
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The latter expression yields dF/dp = p, confirming that p is the density. Clearly, if we 
keep p fixed and take the limit T —* 0, the gauge field becomes singular at the hori- 
zon r = rj\ However, if we eliminate p in favour of p, then the gauge field becomes 
d r Ag = (I^T) 3 V 1 / 2 . At r = tt this is regular for T — > 0. A similar analysis holds 
for the other models analysed in this paper. In contrast, the non-abelian solutions always 
seem to be regular in the T — > limit. 

3.3 Comments on connections with real-world superconductors 

Knowledge of the specific heat in solid state physics is of central importance, as it shows 
what are the relevant degrees of freedom, and how these vary with temperature. In addition, 
this quantity is in a sense "simpler" and more robust than e.g. the resistivity, since, as a 
scalar quantity, it is insensitive to the anisotropic properties of the material. For standard 
metals the specific heat c v behaves as 

c v = cf + $ h = lT + l3T 3 . (3.6) 

The first term, proportional to the Sommerfeld constant 7, originates from an electronic 
contribution and the second term originates from lattice phonons. While the lattice con- 
tribution to the specific heat is unchanged for superconductors, the electronic part varies 
significantly, depending on the type of superconductor (see e.g. [16]). 

The weakly coupled superconductors are successfully described with BCS theory, which 
predicts that the specific heat is exponentially suppressed for temperatures T < T c , i.e. 

= Ae~ B % 

Here 7 is the Sommerfeld constant in the normal phase, and A and B are constants that 
depend on the temperature interval. Therefore, for sufficiently low temperatures the graph 
kx(C/"fT c ) against T c /T is a straight line, which is also experimentally observed (see the left 
panel of figure ^) . The exponential suppression of the specific heat is a natural consequence 
of the presence of a gap, since when T < T c /2, the energy necessary to break a Cooper 
pair is ~ 2A(0), and the number of broken pairs is ~ e - 2A (°)/ fe s' r . Another important 
feature of low T c superconductors is that the specific heat is discontinuous at T = T c since 
at this temperature it is no longer necessary to put energy into breaking of Cooper pairs. 
BCS theory predicts a universal number for this jump, Ac(T c )/7T c = 1.43, which is indeed 
observed in a large number of low-T c superconductors. 

Another class of superconductors are the strongly coupled superconductors which in- 
clude for example high-T c superconductors and heavy fermion compounds. Most of these 
materials are experimentally found to involve higher-spin electric carriers. From this per- 
spective one might expect that these materials should be more similar to non-abelian (spin 
one) condensates considered in this paper. Heavy fermion compounds are very different 
from BCS-like materials when cooled down to the superconducting phase. They have a 
specific heat below T c which is well- approximated with a power law behaviour (see e.g. [19]), 

heavy fermions at T < T c : c v = 70T + kT" , (3.8) 



(3.7) 
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Figure 6: Specific heat of the low-T c BCS superconductor gallium (left) and the heavy- 
fcrmion compound UBei3 (right) in their superconducting phases. Figures taken from [17] 
and [18] respectively. 

where n ~ 2 — 3, and 70 and k are constants. In other words, due to the unusual properties 
of the gap, the specific heat is no longer exponentially suppressed when T < T c . The 
specific heat jump at T = T c , which is observed for the weakly coupled materials, does 
exist for the strongly coupled superconductors as well, although the normalised value of 
the jump is, for many systems, much larger than the value 1.43 predicted by BCS theory 
(see the right panel of figure |6|) . 

When heated above T = T c , a large number of heavy fermion compounds behave as a 
(heavy) fermion liquid, i.e. with specific heat c/T ~ j(T) = const. . However, some of them 
(for example UBei3) have an unusual behaviour and their specific heat is well-described by 



heavy fermions at T > T c : c v ~ — In ( — I or c v ~ T (A « 0.7 — 0.8) . (3.9) 



Keeping all this in mind we can try to compare some of our findings with real super- 
conductors. Firstly, we see that in the superconducting phase, colourful holographic super- 
conductors have power-law behaviour of the specific heat, which is closer to the strongly 
coupled superconductors, not the BCS superconductors. Secondly, all systems exhibit a 
jump in the specific heat, in agreement with real- world superconductors. 

Finally, when T > T c , in all cases the specific heat increases with temperature ac- 
cording to a power law. The precise value of the positive exponents depends on details of 
the system, on whether we work at fixed chemical potential or at fixed density (see the 
previous section) , and finally it depends on whether we work in the Yang-Mills approxima- 
tion (see table j3lj| ) or with the full DBI action [9]. From the experimental perspective the 
specific heat at fixed p is a more interesting quantity, and as we have already explained the 




- 12 - 



Yang-Mills approximation is in this case not appropriate. To analyse the T > T c temper- 
ature range, we thus use the results of [9]. All cases considered in [9] are conformal, such 
that the non-sphere part of the D-brane worldvolume is AdS p +2- They find that c ~ T 2p . 
Hence none of the conformal cases seem to fit the metallic behaviour of heavy fermion com- 
pounds. The specific heat at fixed \x seems more close to real-world superconductors, since 
in the Yang- Mills approximation the D2/D6 system has exponent 1/3, resembling ( p.£| ) 
and D3/D7 has metallic behaviour, while in the DBI approximation the D4/D8 system 
exhibits metallic behaviour [8]. 

In summary, from the perspective of the specific heat, most of the colourful holographic 
superconductors show some qualitative similarities with strongly coupled superconductors. 
In particular, modulo subtleties discussed above, the D3/D7 system looks very similar to 
e.g. the UBei3 heavy fermion superconductor. 



4. Electromagnetic properties 

In this section we turn to a study of the response functions for the superconducting con- 
densates that were constructed in the previous sections. In particular, we focus on the 
computations of AC and DC conductivities. At the end we compare our holographic find- 
ings with qualitative features of real- world superconductors. 

4.1 AC conductivity 

Once a condensate has been found, response functions in this background can be obtained 
by analysing the Green function of the gauge field fluctuations. For the condensate (|2.7| ) . 
one can in principle fluctuate in the direction parallel to the condensate or orthogonal to 
it [4]. Here we choose to consider only fluctuations orthogonal to the condensate, i.e. we 
look at the fluctuation modes in the ifi = A% field. These fluctuation modes satisfy 



Or 



- -d>~rr *22o , a I ~ -<A-33-22/ / |( 1 h2 , / a -<A~00~22 2 , n I A i 

ge fg g d r ip -A^-ge v g g (A\')^-^-ge v g g u ip = . (4.1 



In particular, the Green function is obtained from the boundary data using 

G R {u) = - lim V^e^rt"-^ • (4-2) 
r— >0 

and the AC conductivity then follows from 

a{uj) = . (4.3) 

The fluctuation equation (4.1) exhibits properties just like the condensate equations. 
Following an analysis similar to that for condensate equations (see the appendix), we put 
it in the form 

df(h 3 (f)dri?j - [4(4 1} ) 2 M0 +co 2 h 2 (f)\T- 2 ^ = 0. (4.4) 

From this we conclude that the conductivity a will be a function of the two dimensionless 
parameters jl and p(p) which appear in ( |A.2| ) , as well as the combination uj/T. 
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Figure 7: Real and imaginary parts of the conductivity for the four systems (from 
top to bottom: AdS.4, D2/D6, D3/D7 and D4/D8), plotted for various values of the 
dimensionless chemical potential fi/T (corresponding to the dots in figure)^). 
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In the D3/D7 case, computing the conductivity requires some extra care because of 
the presence of logarithmic terms in the asymptotic expansion of the gauge field fluctua- 
tions [20, 21]. A perturbative analysis yields 

tp = a (l- \u 2 r 2 log(Ar) + ...)+ a 2 r 2 + . . . . (4.5) 

When evaluating the conductivity ( f4.3[ ) this leads to a logarithmically divergent term as well 
as a finite contribution. The logarithmically divergent term can be removed by suitable 
addition of a holographic counterterm. 

We impose ingoing boundary condition on the fluctuation tp at the horizon, which 
means 

i> = (r T - r)~ iau) (l + ai(r T - r) + a 2 (r T - r) 2 + . . .) near r = r T , (4.6) 

where a is a model-dependent real coefficient, \rT, \r^( 2 , \rx, and \y/rr for the AdS4, 
D2/D6, D3/D7 and D4/D8 models respectively. After shooting to the boundary at r = 
we then compute the retarded Green function from ( ^4.2[) . We have analysed the real and 
imaginary parts of the conductivity for all four systems; the results are displayed in figure [7|. 

There is a number of interesting features in both the three and four-dimensional cases. 
First, in comparison to holographic superconductors obtained from AdS4 black holes in 
four-dimensional Abelian Higgs models, such as analysed in [21], we see that for all other 
cases, in the absence of a condensate, the real part of the conductivity does not approach a 
constant. This is in particular true for the two-dimensional, non-conformal D2/D6 system. 
Second, we see that the value of the chemical potential where a gap in the real part of the 
conductivity appears, does not coincide with the critical value of [i where the condensate 
forms. The gap forms at a larger value of and becomes more pronounced as /x increases. 

In all cases we have analysed, there is a strong evidence for the scaling of the energy 
gap with the critical temperature, 

IV 

— constant. (4.7) 
T c 

This is most easily seen from the plots of the conductivity versus the dimensionless ra- 
tio uj/n; an example is given in figure || for the D4D8 system. Following [21], we define uj g 
as the frequency for which Imcr has a minimum (following [21]). The relation between 
the critical temperature T c and the chemical potential ^ is read off from the plots of the 
condensate versus T j fj, in figure ^: the intersection of these curves with the horizontal axis 
yields the value of T/fj, for which the condensate first appears, and hence gives T c for a 



4 An analysis of the non-abelian condensate in the D3/D7 system has previously appeared in [22] but 
that paper ignores the effect of the logarithmic terms in the fluctuation, as witnessed by e.g. the different 
imaginary part of the conductivity. 
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Re [a] 




0.5 1.0 1.5 2.0 -100 L 



Figure 8: The real and imaginary part of the conductivity of the D4/D8 system, ver- 
sus the dimensionless ratio uj/fi, for various values of \x)T (corresponding to the dots 
in figure ||). If we define uj g to be the frequency for which the imaginary part of the 
conductivity has its minimum, this plot shows that oj g ~ fi. The other systems show 
qualitatively similar behaviour. 



given /j,. 5 In this way we obtain 



AdS 4 : 




s 0.13, 


LOg/T c P 


d 8.4 


D2/D6 : 




s 0.14, 


Wg/Tc P 


t 8.0 


D3/D7 : 


T c /fi p 


s 0.16, 


w 9 /T c p 


i 8.1 


D4/D8 : 


T c /^p 


a 0.21, 


w g /T c p 


d 6.7 



(4- 



In the conformal cases analysed thus far this ratio turned out to be around 8 to within a 
few percent. This is in contrast to the BCS values which are anywhere between 3.5 at weak 
coupling and 4 in the naive strong coupling limit (see e.g. [23]). The non-conformal cases 
give results which are higher than the BCS result too, with the D4/D8 value deviating 
substantially from 8. Our results for the non-abelian condensate in the AdS 4 black hole 
background is similar to the one for the AdS 4 Abelian Higgs model of [21]. 

One final way to compare the holographic superconductors with experimental data on 
superconductors is to express the results in terms of the real resistivity, defined as 

R = (Rea) 2 + (Ima) 2 ' (4 ' 9) 

Although plots of R versus uj for the AC behaviour are not very common in the literature, 
we have provided these for completeness in figures ||. We will have more comments on the 
behaviour of the DC resistivity in the next section. 

4.2 DC conductivity 

It has recently been noted that the u = limit of the real part of the conductivity can be 
expressed in terms of quantities at the horizon, rather than at infinity [24]. By expressing 



5 In reading off the numerical coefficient, it should be noted that the variable along the horizontal 
axes of these plots is qT/fi, with the proportionality constant q chosen such that this reduces to 
when tt = L = 1. 
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Figure 9: Resistivity as a function of frequency (left to right, top to bottom: AdSzt, 
D2/D6, D3/D7 and D4/D8), for various values of /i/T (corresponding to the dots in 
figure |||). Lighter curves correspond to larger chemical potential. 



the conductivity in terms of the canonical momentum, and noting that it satisfies a trivial 
flow equation in the radial direction. The conductivity can be expressed in terms of a 
quantity evaluated at the horizon. A generalisation of their argument to include a dilaton 
reads 



Reo-DC °c \ j — —g xx e 

9rr9tt 



r=r t 



(4.10) 



(where the proportionality involves an overall coupling constant dependence which will not 
be important here). For the cases analysed here, this expression yields Re a oc i, rvs, T 
and T 2 respectively 6 . It is important to note that, except for the conformal AdS4 case, these 
are all increasing functions of temperature, quite in contrast to e.g. ordinary metals. We 
should also note that while the result of [24] was derived in the absence of any condensate, 
it is also relevant for the abelian condensate, since in the Yang-Mills approximation the 
fluctuation equations decouple from the abelian condensate. 

This behaviour of the DC conductivity at zero condensate is indeed easily verified to 
hold true for our numerical results in the uj — > limit, thus providing a quick test of our 
analysis in the abelian case. In fact, one can again use the properties of equation ( }4.4| ) 
explained in the appendix to analyse the behaviour of a as a function of T. First, one 
notes that in r coordinates the Green function (|4.2f) is given by an overall power of TL 



The D3/D7 behaviour agrees with [25] and is the same as found for the Abelian Higgs model in [21]. 
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Figure 10: The real and imaginary part of the DC conductivity (or rather, the AC 
conductivity at fixed small frequency) as a function of temperature, and the resistivity, 
for the D3/D7 system. For temperatures T > T c the condensate becomes abelian and 
we recover the Re a oc T result of [24, 25], which is the large-temperature limit of [20]. 
For T < T c there is a regime for which Rec is not yet zero, but shows a polynomial 
dependence on the temperature. 



multiplied with a function of uj/T. These factors of TL are 



AdS 4 : 


G R (u 


T) 


= (TL)xG R (u/T, fi /T), 


D2/D6 : 


Gr(oj 


T) 


= (TL)I 


xGr(lo/T,h/T) 


D3/D7 : 


G R (u 


T) 


= (TL) 2 


x G R {u/T,n/T) 


D4/D8 : 


G R (u 


T) 


= (TL) 3 


x G R {u/T,n/T) 



(4.11) 



The only bare oj dependence enters through the denominator of ( |4.3| ) . 

Once the non-abelian condensate is turned on, the analysis of [24] can no longer be 
applied. The coupling between the fluctuation and the gauge field component A x implies 
that it is now no longer true that the radial evolution of the canonical momentum is trivial. 
Hence it is no longer possible to express the conductivity in terms of a quantity evaluated 
at the horizon. Nevertheless, it is possible to analyse the temperature dependence of the 
conductivity in the u — > limit numerically. The scaling properties imply that ji/T oc uj/T. 



Figure |10| displays the behaviour of the conductivity and resistivity in this non-abelian 
phase. 7 

Generically, we see that superconductors of the holographic type exhibit a behaviour of 
the resistivity which is similar to that of semi-conductors, i.e. R decreases with temperature 
above the critical temperature. A metallic behaviour has so far only been found for small 
temperatures in the model of [20]. 

4.3 Comments on connections with real-world superconductors 



We have already seen in section that weakly coupled and strongly coupled superconduc- 
tors exhibit very different behaviour for their specific heat. As far as the electric resistivity 
of these materials is concerned, while all of them have zero resistivity for T < T c , when 



7 The metallic Recr oc 1/T 2 behaviour at small temperatures, found in [20], is invisible in the Yang- 
Mills truncation. It is, however, unclear whether it will show up before the transition to the non-abelian 
condensate takes place. 



- 18 - 




Figure 11: Resistivity of the heavy fermion compound UBei3. Figure taken from [26]. 



T > T c , their behaviour again differs significantly. The weakly coupled superconductors 
mainly behave as "metals" above T c (i.e. with a resistivity that grows with T as a poly- 
nomial of T). On the other hand, just above T c strongly coupled superconductors behave 
either as metals or as semiconductors, with a resistivity that decreases with increasing 



temperature (see figure 11). 

Just like real superconductors, all holographic models exhibit zero conductivity when 
T < T c . Interestingly, the temperature at which the resistivity becomes strictly zero is 
lower than the temperature T c at which the non-abelian condensate becomes the true 



ground state. This is most clearly visible in figure 10. The small temperature range 
To < T < T c corresponds to a transition region where the resistivity very quickly drops 
from the normal-phase value to zero (the drop is more clearly visible in the real part of the 
conductivity, as the resistivity is suppressed by the large imaginary part). 

When T > T c only Aq 7^ 0, and hence only the abelian condensate p determines the 
behaviour of the system. The Yang-Mills truncation leads to a conductivity which is inde- 
pendent of the condensate (or chemical potential), due to the fact that the electromagnetic 
fluctuations fully decouple from the condensate in the abelian approximation. In this ap- 
proximation, the system behaves as a semi-conductor, with a resistivity which decreases 
with increasing temperature (see figure [Toj ) . The DBI approximation, on the other hand, 
couples the fluctuations to the condensate, and leads to a condensate-dependent result for 
the conductivity [20], 

* 2 = aT 2 + ^, (4.12) 

where a and b are numerical constants and A is the 't Hooft coupling. The first, condensate- 
independent term originates from charge carriers that are thermally produced in charge 
neutral pairs, while the second term originates from the charge carriers introduced by the 
abelian condensate. The first term was visible in the Yang-Mills approximation and leads to 
semiconducting behaviour (an increase of temperature enhances thermal pair production), 
while the second term describes metallic behaviour. 



- 19 - 



The system thus behaves as a metal below the temperature T < T*, 

r.-(S)". 

which is governed by a' corrections to the Yang-Mills action. As the system is cooled, the 
question is whether metallic behaviour will kick in before the superconducting phase or 
not. Obviously, this depends on whether T c is smaller or larger than T*. The requirement 
that > T c , i.e. that the system behaves as a metal above T c , will put constraints on 
A, and determine if it is compatible with the requirement that A <C 1. Of course, one 
should keep in mind that our analysis for T < T c has been done in the (non-abelian) Yang- 
Mills approximation, while the analysis for T > T c was done in the DBI approximation. 
Therefore, the first step one should take is to see how higher-derivative corrections could 
change the critical temperature. We leave this issue for future investigation. 
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A. Appendix: computing temperature dependence 



While the coupled system of equations ( |2.9| ) is complicated and does not generically admit 
analytic solutions, it does have a useful scaling symmetry (both in the conformal and 
nonconformal cases), which simplifies our analysis. When considering a non-conformal 
systems we deal with three dimensionful parameters, the scale L, the temperature T and 
the chemical potential fi. One might thus expect that the dimensionless condensate is a 
function of the two dimensionless combinations TL and \xL. However, the equations ( |2.S| ) 
allow us to reduce this dependence to only one dimensionless combination fi/T, with all TL 
dependence following from general arguments. This is a consequence of a scaling symmetry 



of (2^9), even when the metric (2.2) does not have conformal symmetry. 

To see this, let us introduce a dimensionless radial coordinate r = r/r^. In this 
coordinate the equations that determine the condensate read 

d f (/ii(f) fyA®) = 4(4 1) ) 2 4 3) ^(f) T' 2 , 

(A.l) 

d f (/i 3 (f) ^4°) = 4(4 3) ) 2 4 1} Mf) T- 2 , 

where the hi(r) are functions of f only. That is, all T-dependence sits in the explicit T~ 2 
factor of the terms on the right-hand side, and there is no bare dependence on L anymore. 
Solutions to the equations now take the form 

4 3) = n + p r a + . . . = tU + pr a + . . . ) =: TAf , (A.2) 
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and similar for A^. (see e.g. (A. 4) for the D4/D8 case). Here a is a model-dependent 
constant and we have defined 



p := 



— and 
T 



P 



PVT) 



(A.3) 



terms of the Aq and A^' variables, we can then conclude that the dimensionless quantity 
p only depends on p, not on the dimensionless combination TL. For the D4/D8 system, 
for instance, the expansion reads 



f(i) 



From the fact that the equations (|A.1|) are independent of T and L when expressed in 



D4/D8 : A { q ] 



P 



T T 5 / 2 (TL) 3 / 2 



f 3/2 + 



P 



pr 3 / 2 + 



(A.4) 



So we see that even when the theory is not conformally invariant 8 , the scaling symmetry 
of the equations allows us to determine the dependence of the dimensionless condensate 
on TL. We also note that if one were to include higher derivative terms into the equations of 
motion (which for example follow from the DBI action), then the above symmetry would 
not necessarily continue to be present. In practice, we have used the scaling symmetry 
to verify numerical stability of the condensate solutions: plots of p against p should be 
independent of TL. 

The scaling argument also allows us to determine e.g. the temperature dependence 
of the action by changing p instead of rj-. The idea here is to compute the action in r 
coordinates at a fixed value = 1. This will yield S(p, T) = r^S(p/T, tt = 1) for some /?, 
with the overall power of tt simply arising from the coordinate transformation r = r^r. 
We will see this at work in section |3.lL 



s In the conformal D3/D7 case, this story simplifies and the dimensionless condensate p/T is independent 
of TL, but the argument presented here shows that the story is somewhat more complicated in the non- 
conformal cases, where the naive dimensionless combination p/T' Q+1 - t still depends on TL, albeit only in 
the simple way dictated by (|A.3|). 
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